Searching for statistical equilibrium in a dynamical multifragmentation path 
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A method for identifying statistical equilibrium stages in dynamical multifragmentation paths as 
provided by transport models, already successfully tested for the reaction ^^^Xe+^^^Sn at 32 MeV/u 
is applied here to a higher energy reaction, ^^®Xe+^^^Sn at 50 MeV/u. The method evaluates 
equilibrium from the point of view of the microcanonical multifragmentation model (MMM) and 
reactions are simulated by means of the stochastic mean field model (SMF). A unique solution, 
corresponding to the maximum population of the system phase space, was identified suggesting that 
a huge part of the available phase space is occupied even in the case of the 50 MeV/u reaction, 
in presence of a considerable amount of radial collective flow. The specific equilibration time and 
volume are identified and differences between the two systems are discussed. 

PACS numbers: 24.10.Pa; 25.70.Pq 
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I. INTRODUCTION 

From more than 20 years, statistical equilibrium was a 
basic ingredient in many theoretical [l|, d, H, d, [1| and 
experimental nuclear multifragmentation studies. This 
was motivated by the good agreements obtained over the 
time between various observables related to the asymp- 
totically resulted fragments and various statistical mul- 
tifragmentation models. 

However, the agreement in the asymptotic stage be- 
tween experiments and statistical models predictions is 
only a necessary condition for equilibrium. For fully 
"proving" equilibrium one has to have direct access to 
the primary decay stage of the reaction. In fact, differ- 
ent primary configurations could lead to the same final 
results, because of compensative effects between primary 
and secondary emission mechanisms. Experimentally it 
is difficult to "measure" the primary decay stage: this 
signal is distorted by the secondary particle emission. 
An elegant solution is to perform a statistical analysis 
on the primary decay stage as predicted by a dynamical 
model yielding results in good agreement with experi- 
mental data. One of such models is the Stochastic Mean 
Field (SMF)[i,0,i. 

In Ref. a statistical analysis was performed on the 
dynamical path provided by SMF with the microcanon- 
ical multifragmentation model (MMM) 5]. There, an 
effective method for identifying the equilibrated stage in 
the dynamical path was proposed and successfully tested 
for the i29Xe-H"9Sn at 32 MeV/u reaction. A fully equi- 
librated source was identified at 140 fm/c, with a corre- 
sponding volume of 3.4 Vq (where Vq is the volume of the 
source at normal nuclear density). Herein we apply the 
same analysis to the higher energy reaction ^^^Xe-|-^^^Sn 
at 50 MeV/u as described by the SMF model. In this 
way one can discuss the possible occurrence of statisti- 
cal equilibrium even for a more explosive system, having 
a considerable amount of radial collective fiow. Impor- 
tant differences between the two reactions are supposed 
to appear. Since at 50 MeV/u the emitted fragments are 



generally smaller, one can find a large variety of source 
configurations that can fit the fragment partitions. How- 
ever, we will show that the analysis of the corresponding 
fragment kinetic properties allows to clearly identify a 
unique solution. 



II. BRIEF REVIEW OF THE EMPLOYED 
MODELS 

Let us first give a brief description of the two models 
involved in this analysis. According to the SMF theory 
J5J, during the expansion phase that follows the initial 
coUisional shock the system encounters volume (spinodal) 
and/or surface instabilities, that lead to its break-up into 
many pieces. As in mean-field approaches, the system 
is represented by the one-body distribution function in 
phase space, that evolves according to the self-consistent 
nuclear (+ Coulomb) potential and to the stochastic col- 
lision integral, that accounts for the residual two-body 
interaction. In the approach considered, fluctuations are 
introduced in an approximate way Q . Their amplitude is 
essentially determined by the degree of thermal agitation 
present in the system. Then fluctuations are amplified 
by the unstable mean-field, leading to the formation of 
fragments, whose properties reflect the structure of the 
most unstable collective modes. Several multipoles are 
excited with close probabilities, leading to a large variety 
of fragment configurations. ^lOj . 

Concerning more technical aspects, such as fragment 
recognition, we follow a coalescence procedure of the one- 
body density. Fragment excitation energies are calcu- 
lated adopting the local density approximation, by sub- 
tracting the Fermi motion (associated with the local den- 
sity) from the fragment kinetic energy (taken in the frag- 
ment reference frame) [13, [HI ■ 

The MMM model [a| describes the break-up of a statis- 
tically equilibrated nuclear source defined by the param- 
eters: mass number (A), atomic number (Z), excitation 
energy (E) and freeze-out volume {V). The model as- 
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sumes equal probability between all configurations C : 
{Ai, Zi, ei,ri,Pi, i = l,...,iV} (the mass number, the 
atomic number, the excitation energy, the position and 
the momentum of each fragment i of the configuration 
C, composed of N fragments) subject to microcanoni- 
cal constraints: J2i = A J2i ^» = J2i P» = 0- 
ri X Pi = 0, E - constant. The level density of a 
given fragment (entering the statistical weight of a con- 
figuration) is taken to be of Fermi-gas type adjusted with 
the cut-off factor exp(— e/r): p(e) = po{e) exp(— e/r) [l^l 
that counts for the dramatic decrease of the lifetime of 
fragment excited states respective to the freeze-out spe- 
cific time as the excitation energy increases. MMM can 
work within two freeze-out hypotheses: (1) fragments are 
treated as hard spheres placed into a spherical freeze-out 
recipient and are not allowed to overlap each-other or 
the recipient wall; (2) fragments may be deformed and 
a corresponding free-volume expression is approaching 
the integration over fragment positions ^12] . Hypothe- 
sis (2) is more adequate for the present study since it 
allows exploration of higher density configurations. In- 
deed, our dynamical studies seem to indicate that a good 
amount of statistical equilibrium is reached at a pre- 
fragment level, i.e. when the produced fragments are 
not fully separated yet [§]. This condition cannot be 
always matched by the space configurations related to 
the hypothesis (1), which presents severe constraints for 
low frcczc-out volumes. Moreover, as resulting from p^ . 
at relatively small freeze-out densities the results of the 
two approaches roughly coincide. This is why we per- 
fom the present investigation within hypothesis (2). The 
model employs a Metropolis-type simulation for estimat- 
ing the average value of any system observable (see Refs. 
y for more details). While MMM includes a secondary 
decay stage, this stage is not necessary for the present 
analysis because we aim at a direct investigation of the 
(primary) break-up stage. A feature of this model, partic- 
ularly important for the present study is the possibility of 
including flow in the primary decay stage. We include a 
standard flow velocity profile v = wo(r/i?)" where r 
is the distance of a given fragment from the system cm. 
and R is the freeze-out recipient radius. In MMM the 
total flow energy is microcanonically conserved. There- 
fore, Wo is evaluted for each fragmentation event from the 
condition of conservation of the system flow energy. The 
flow exponent a will be varied between the values 1 and 
2. 



III. STATISTICAL ANALYSIS 

Using MMM we perform the analysis proposed in Rcf. 
in order to identify a possible statistical equilibration 
stage in the dynamical path of the reaction ^^^Xe+^^^Sn 
at 50 MeV/u as provided by SMF. Within SMF the dy- 
namical freeze-out time is defined as the time when the 
fragment formation process is over and therefore (inter- 
mediate mass) fragment multiplicities do not evolve any 



longer. This time was identified to be around 200 fm/c. 

In order to wash-up pre-equilibrium effects, which 
should appear in the dynamical simulation, only interme- 
diate mass fragments (IMF) (i.e. fragments with Z > 3) 
are considered from both SMF and MMM simulations. 
Due to the large Coulomb repulsion among primary frag- 
ments and the possible presence of a collective flow, it is 
reasonable to assume that from the equilibrated freeze- 
out, i.e. the time when the system could have reached 
full phase space occupation, to the readily identified dy- 
namical freeze-out the major difference will be the volume 
in which fragments are located. In other words we con- 
sider that all source variables except volume are roughly 
preserved. The fragment properties included in our anal- 
ysis are connected to relevant inclusive observables that 
can be predicted by both models: fragment multiplicity 
and distributions and fragment internal excitation en- 
ergy. We will fit these properties, but not the volume. 
For finding the best fit we have to minimize the follow- 
ing error function [^: 

£l = {3[fi{Abound))+f{{Zbound))] 



+ 



J2 f[{dN/dNiMF)]/ J2 1 



.NI^ 



+ fi{eIMF))+J2fi{Z.r...^))m/9 (1) 
1=1 

where (•) stands for average. Abound and Zbound are 
the bound mass and charge (sum of the mass number 
and, respectively, atomic number of all IMF's from a 
given event), Njmf is the number of IMF's, eiMp is 
the fragment excitation energy per nucleon and Zmaxi 
with i = 1,2,3 are the largest, second largest and third 
largest charge from one fragmentation event. Further, 
f{x) = \2{xs — Xd)/{xs -f Xd)\, where the indexes s and 
d stand for "statistic" and "dynamic" . 

For finding the coordinates of the minimum of £i we 
variate the MMM parameters A, E, V and r in wide 
ranges thus constructing a four-dimensional grid. The 
ranges are A : [170,210], E : [4,8.5] MeV/u, V/Vq : 
[2, 14], T = 20, 50, 100, oo MeV. The source is considered 
to have the A/Z ratio of the ^^^Xe-|-^^^Sn reaction. Like 
for the case of the 32 MeV/u incident energy for a given 
value of T one gets a unique minimum for £i. This is 
illustrated in Fig. 1 where cuts in £i corresponding to 
Amin and i^inin are represented for two values of r: 50, 
oo MeV (where the index 'min' refers to the coordinates 
of the minimum value of However, apart from the 
^■^^Xe-l-^^^Sn at 32 MeV/u reaction, in the present case 
one cannot get an unambiguous source identification only 
by minimizing £i. Indeed, it is possible to find a solution 
for any value of the considered t and, as resulting from 
Fig. 2, the curve representing the minimum value of £i 
versus T is rather flat. In other words, primary fragment 
distributions and their internal excitation energies are 
not sufficient for proving the occurrence of equilibrium. 

So, in order to better identify the equilibrated source we 
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FIG. 1: Contour plots of the error function £i [see eq. (1)] for 
r = 50 MeV (left column) and t — oo MeV (right column). 
Cuts in the (V/Vo,E) (upper panel) and in the {A,V/Vo) 
(lower panel) planes corresponding to the minimum £i 
coordinates. Darker regions correspond to smaller £i; units 
are relative. 

have to analyze the reaction kinematics too. More specif- 
ically, we fit the fragment kinetic distributions and the 
positions of the three highest charged fragments corre- 
sponding to the dynamical freeze-out. The error function 
employed in this fit is: 
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5]/((M^)))/36 + ^/((r(Z.„ax.)))/3 



.2=5 



/2 (2) 



{k denotes the fragment kinetic energy and r denotes its 
position; in practice, for increasing the accuracy of the 
£2 estimator, a smoothed k{z) curve was used and the z 
interval was spanned with a step of 5 units). But how 
can one get "statistical" information about the dynam- 
ical freeze-out in order to perform this fit? One simply 
has to propagate the primary fragments in their mutual 
Coulomb field, taking into account also the presence of 
collective flow effects, from their freeze-out positions as 
generated by MMM up to Vjmf = 16.73 Vq which is 
the volume corresponding to the dynamical freeze-out. 
(Herein we use the notations from Ref. [9]: V is the 
volume of the smallest sphere which totally includes all 
fragments; V the volume of the smallest sphere which 
totally includes all fragments and has the center located 
in the center of mass of the system) . Such propagations 
were performed for sources corresponding to the mini- 
mum £1 coordinates for each of the considered values of 
r. The flow energy, Efi = ^nnvf {mi and Vi being 
respectively the mass and the flow energy of the frag- 
ment i), and the flow exponent were varied in the ranges: 
Efi : [0, 1.6] MeV/u and a/; : [1, 2]. An unique minimum 
of the £2 function was found for each of the considered 
values of r. One can observe this behavior in Fig. 3 
for the case r = 00. The obtained dependence of the £2 
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FIG. 2: The error functions E\, £2 and £=£±+£2 versus 1/r. 

minima of r is given in Fig. 2. This time we observe a 
monotonic dependence of f2min of The global error 
function £ = £1 -\- £2 exhibits a well defined minimum as 
well (from now on we denote by £ this minimum value). 
This minimum corresponds to t = c». 

It should be noticed that the full set of inclusive ob- 
servables considered in the total estimator, £, is quite 
exhaustive. Most of them coincide with the observables 
usually exploited in the comparison with experimental 
data. 

For the sake of completness, we test the £2 estimator 
in the case of the 32 MeV/u reaction, which was already 
analyzed in Ref. [9j. There it was shown that fragment 
kinetic energies and positions were compatible with the 
source identified only with £1 for a zero flow energy. In 
agreement with Ref. [i^ results, the minimum of the £2 
versus r curve sharply points the value r = 00, as £1 does, 
see Fig. 4. In the {Efi, ay;), t — 00 plane £2 minimizes 
inside the Efi—0 MeV region, as shown in Fig. 5, which 
confirms the findings of Ref. ^9J . The sharper minimum 
of the £ = £1 + £2 estimator observed for the 32 MeV/u 
reaction compared to the 50 MeV/u situation indicates a 
more advanced degree of statistical equilibrium reached 
in the 32 MeV/u reaction. 

Finally, it is worth noticing that since hypothesis 1) 
of the MMM model induces artificial geometrical con- 
straints at small volumes, it is easy to infer that, for 
example, in the case of 32 MeV/u one would not find 
any physically meaningfull minimum (i.e deep enough to 
insure a good reproduction of the data under study) . On 
the other hand, for the 50 MeV/u reaction, the system 
being less affected by these constraints, the result would 
be closer to the evaluations obtained within hypothesis 
2). 



IV. INTERPRETING THE RESULTS 

Some discussion is now in order for the 50 MeV/u re- 
action case: while mathematically speaking the absolute 
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TABLE I: MMM: r = 50 MeV, oo versus SMF results for 
various system observables. (eiMp) is given in MeV/u units. 



0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 

Efi (MeV/u) 

FIG. 3: Contour plot of the function £2 in the plane {Efi, Ufi) 
corresponding to r = 00. Units are relative; darker regions 
correspond to smaller values of £2- 
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FIG. 4: The error functions £i, £2 and £—£i+£2 versus 1/r 
for the case of the reaction Xe+Sn at 32 MeV/u. 
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FIG. 5: Contour plot of the function £2 in the plane {Efi, Ufi) 
for the case of the reaction Xe+Sn at 32 MeV/u. Units 
are relative; darker regions correspond to smaller values of £2- 
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minimum corresponds indeed to r = 00, one may ob- 
serve a plateau-like behavior of this quantity starting at 
T = 50 MeV. It means that all values of r between 50 
MeV and 00 are almost equally good solutions. An illus- 
tration of this fact is given in Fig. 6 where charge, num- 
ber of intermediate mass fragments, kinetic energy and 
three highest charged fragment position distributions cor- 
responding to T =50 MeV and 00 calculated with MMM 
are compared with the SMF results. The correspond- 
ing source parameters are: (1) r = 50 MeV, A=183, 
Z=76, E^6.6 MeV/u, F/Vo=5.6, F/Vb=6.79, Efi^O.6 
MeV/u, a/i=1.2, ip=58.8 fm/c; (2) r = 00, A=183, 
Z=76, E=6A MeV/u, V/Vo=9.6, V/Vo^U.78, Efi^l 



MeV/u, a/; = 1.4, tp~lA.8 fm/c. (Here tp denotes the 
propagation time, i.e. the time needed to reach, starting 
from the configuration predicted by MMM, the dynam- 
ical freeze-out configuration). The resemblance between 
the two cases transpares also from Table I, were the av- 
erage value of some observables involved in the construc- 
tion of £, evaluated by MMM for r = 50 MeV, 00 are 
compared with the corresponding SMF ones. In terms 
of time along the dynamical evolution, the r = 50 MeV 
and 00 stages correspond respectively to approximately 
140 fm/c and 185 fm/c. Note that for the 50 MeV/u 
reaction case, the agreement betwen dynamical and sta- 
tistical radial distributions of the three largest fragments 
is weaker compared to the case of the 32 MeV/u reaction 
- see Fig. 7 (Ref. Q). 

But how to interpret these results? Which is the real 
equilibration time? Turning back to Fig. 2 we ob- 
serve that the major change in £ versus r occurred at 
r = 50 MeV. An interpretation would be that while a 
major part of the system phase space has already been 
spanned at t =140 fm/c (i.e. r = 50 MeV), slight ad- 
justments towards equilibrium are further achieved until 
t =185 fm/c (i.e. t = 00). 

Two facts are strongly supporting this hypothesis: 
(1) As done for the 32 MeV case [9)], one can per- 
form, along the dynamical evolution, an analysis in terms 
of pre- fragments, i.e. try to identify density bumps 
(whose density is larger than the average) that will fi- 
nally develop into the actual fra gme nts observed at the 
dynamical freeze — out time ll|, 1^ . The SMF number 
of pre-fragments saturates at 140 fm/c (see Fig. 8). This 
is in agreement with Ref. 9| findings where the equilibra- 
tion time for the ^^^Xe-f^^Sn at 32 MeV/u was found to 
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FIG. 6; From top to bottom; charge, number of intermediate 
mass fragments (IMF), kinetic, position of the first three 
highest charged fragment distributions (the position distri- 
butions are normalized by 4-Kr^). Left column: r = 50 MeV 
MeV; right column: t — oo MeV. MMM results (squares - 
first three rows from the top; solid line - last row) are plotted 
in comparison with the SMF ones. Results are corresponding 
to the i^^Xe+^'^Sn at 50 MeV/u reaction. 




r (fm) 

FIG. 7: Distributions of the position of the first three highest 
charged fragments (normalized by Airr^). MMM results 
(solid line) are plotted in comparison with the SMF ones 
(dashed lines). Results correspond to the ^'^^Xe+'^^^Sn at 32 
MeV/u reaction [1]. 
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FIG. 8: (Color online) Dynamical results: Evolution of IMF 
(circles) and pre-fragment (squares+dashes) multiplicities in 
time. 
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FIG. 9: Spatial distribution of test-particles in SMF corre- 
sponding to the moments of time 140 fm/c (left) and 200 
fm/c (right) for an arbitray fragmentation event. 

coincide with the pre-fragment number saturation point. 

(2) While at 140 fm/c the system is still in a pre- 
fragment stage, and therefore there is strong interaction 
between the fragment seeds, at 185 fm/c fragments are 
almost completely formed and separated. This is clearly 
illustrated by Fig. 9 where fragment spatial configura- 
tions yielded by SMF, corresponding to the two moments 
of time are plotted. From this perspective, the further 
adjustments towards equilibrium taking place until 185 
fm/c are achievable through a strongly decreasing nuclear 
inter-fragment interaction and Coulomb repulsion. 

As a test for the validity of the analysis in terms of 
pre-fragments, we represent in Fig. 10 the distribution 
of the first three highest charged (pre) fragments for two 
moments of time: t =140 fm/c and t = 200 fm/c. As 
it can be observed in Fig. 10 there is a perfect agree- 
ment between the distributions corresponding to the con- 
sidered moments of time. This fact proves that the 
(pre)fragments partitions remain practicaly unchanged 
during the expansion from the statistical freeze- out to 
the dynamical freeze- out, nicely supporting the working 
hypothesis of our analysis. 
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FIG. 10: Distributions of the first tliree liigfiest charged 
pre-fragments corresponding to two moments of time of the 
reaction Xe+Sn at 50 MeV/u. Full lines correspond to 
t — 140 fm/c; dashed lines correspond to t = 200 fm/c. 



While fairly acceptable, the general agreement between 
the SMF results and the MMM ones for ^^Sxe+^^Sn at 
50 MeV/u reaction is weaker compared to the one cor- 
responding to the ^^^Xe+^^^Sn at 32 MeV/u reaction. 
In the latter case the very good agreement between sta- 
tistical and dynamical results for all the considered ob- 
servables, is a strong indication of the system full equili- 
bration. In particular, the distribution of the position of 
the three highest charged fragments, that is reported in 
Fig. 7 for the 32 MeV/u case, is well reproduced, while the 
agreement is worse at 50 MeV/u (compare the bottom 
part of Fig. 6 and Fig. 7). Part of this difference is due to 
the approximations involved in both approaches. In fact, 
a 20% higher flow energy is observed at the dynamical 
freeze-out time as compared to the statistical one, sug- 
gesting some topological differences between the freeze- 
out partitions of the two approaches. The other part of 
it may be simply due to the fact that in the case of the 
50 MeV/u reaction the system phase space is only al- 
most spanned. This could be attributed to the fact that 
the present reaction is more violent compared to the 32 
MeV/u one and, therefore, "more" dynamically driven. 
In this respect, see for example the more evident hol- 
low configuration corresponding to the dynamical events 
(Fig.6-bottom panels and Fig. 9). 



V. FINAL REMARKS 

Summarizing, the statistical analysis on the SMF dy- 
namical path, started in Ref. Q, was extended to the 
i29Xe-H"9Sn at 50 MeV/u reaction. While most of the 
system configuration space appear to be already spanned 
at 140 fm/c, slight adjustments towards equilibration are 
further achieved until 185 fm/c. 

The following conclusions about the statistical equili- 
bration mechanism in dynamical paths can be drawn: 1) 
equilibrium occurs while system constituents are still in 
interaction; 2) an external radial flow constraint is neces- 
sary at 50 MeV/A. For both ^^Sxe-hH^Sn at 32 MeV/u 
[H and 50 MeV/u reactions (most) of the equilibrium is 
already reached in a pre-fragment stage of the system. 
For both reactions this stage was found to coincide with 
the saturation of the pre-fragment number. Remarkably, 
this time is identical for both 32 MeV/u and 50 MeV/u 
reactions: phase space is spanned after the system has 
spent approximately 50 fm/c inside the unstable region 
of the liquid-gas phase diagram and corresponds to about 
140 fm/c after the beginning of the reaction. Indeed, this 
time is related to the time scale of the instability growth 
towards the pre- fragment formation (6l.[l0|. However, the 
corresponding freeze-out volume is 5.6Vo, much larger 
than the one corresponding to the 32 MeV/u reaction, 
due to the larger compression-expansion effects present 
in the 50 MeV/u reaction case. Further, it should be no- 
ticed that, while in the 32 MeV/u reaction case the equili- 
brated configuration corresponds to r = oo (see Ref. [9] ) , 
here, for the 50 MeV/u reaction, we find t = 50 MeV. 
The physical meaning of this parameter, and its relation 
to the density and temperature conditions reached along 
the reaction path, need to be further investigated. 

The present analysis reveals a qualitative difference be- 
tween the 32 MeV/u and 50 MeV/u reactions: While in 
the 32 MeV/u case @ fragment partitions and excita- 
tion energies impose strong constraints already identi- 
fying the equilibrated source, in the 50 MeV/u case, a 
subsequent kinematic analysis is necessary. This may be 
due to the higher amount of liquid-like phase in the 32 
MeV/u reaction (see e.g. the more pronounced shoulder- 
like shape in the 32 MeV/u charge distribution of Ref.[9]), 
adding extra fragment-size related constraints to the 32 
MeV/u events compared to the 50 MeV/u ones. This 
detail should be instructive for multifragmentation stud- 
ies dealing with identifying the equilibrated sources by 
fitting the data with statistical models. 

Finally, the present results may be useful for the exper- 
imental nuclear thermodynamics studies where volume is 
a key variable for locating the process in the phase dia- 
gram. 
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